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A model is investigated which displays a picture of the symmetry energy as an energy of rotation 
in isospace of a Cooper pair condensate, briefly "superfluid isorotation" . The Hamiltonian is isobari- 
cally invariant and has a one- and a two-nucleon term, where the two-nucleon interaction is composed 
of an isovector pairing force and an interaction of isospins. It is analyzed in the Hartree-Bogolyubov 
plus Random Phase approximation. The Hartree-Bogolyubov energy minus Lagrangian multiplier 
terms proportional to the number of valence nucleons and the z component of the isospin is shown 
to be locally minimized by a product of neutron and proton Bardeen-Cooper-Schrieffer states. The 
equations of the Random Phase Approximation (RPA) can be reduced to independent equations 
for two-neutron, two-proton, and neutron-proton quasiparticle pairs. In each of these spaces, they 
have a Nambu-Goldstone solution due to the global gauge invariance and isobaric invariance of the 
Hamiltonian. Except for the Nambu-Goldstone solutions, the RPA solutions are independent of the 
strength of the isospin interaction. If, in one space, the pertinent single-nucleon spectrum has a 
particle-hole symmetry, the RPA solutions are twofold degenerate except for the Nambu-Goldstone 
solution and one more solution. In an idealized case of infinitely many equidistant single-nucleon 
levels, the one-nucleon term in the Hamiltonian and the isospin interaction contribute terms in the 
symmetry energy quadratic in the isospin T. The pairing force and the two-neutron and two-proton 
RPA correlation energies do not contribute. The contribution of the neutron-proton correlation 
energy is dominated by the Nambu-Goldstone solution, which gives a linear term that makes the 
total symmetry energy proportional to T{T + 1). The rest of this contribution is negative and can 
be written as the difference of two terms of the form ^y(clT)^+lP — b. Observations reported from 
Skyrme force calculations are discussed in the light of these results. Calculations with deformed 
Woods-Sax;on single-nucleon levels give results similar to those of the idealized case. In calculations 
for the mass numbers j4 = 56 and A — 100 with spherical Woods-Saxon levels, the promotion of 
nucleons across magic gaps in the single-nucleon spectrum and the onset of superfluidity with the 
departure from magicity give rise to large linear terms in the symmetry energy. The calculations 
with Woods-Saxon single-nucleon levels reproduce surprisingly well the empirical symmetry energy. 
An experimental signature of superfluid isorotation is discussed. 



I. INTRODUCTION 

The concept of a symmetry energy originates in 
Weizsacker's early attempt [1] to construct a formula 
for the nuclear binding energy. Guided by Majoranas 
ideas [2] as to the nature of the internucleon force, 
Weizsacker suggests an expression for the binding en- 
ergy of a doubly even nucleus which, except for an elec- 
trostatic term, is symmetric in the numbers N and Z 
of neutrons and protons. If Bs{N,Z) is the symmet- 
, ric part of such an expression, the symmetry energy is 
Ss(yl/2, A/2)-Bs{N, Z) with A = N + Z. As discussed, 
for example, by Bohr and Mottelson [3], it carries infor- 
mation on basic aspects of the nuclear structure. Models 
of the nucleosynthesis in core-collapse supernovae rely 
on estimates of the masses of nuclei inaccessible to ex- 
periment. The accuracy of such estimates depends on a 
valid understanding of the origin of the symmetry energy, 
and data on the abundancies of nuclides in the nearby 
universe may in turn constrain nuclear models, as dis- 
cussed in recent reviews by Arnould and Goriely [4] , and 
Arnould, Goriely, and Takahashi [5]. In infinite nuclear 
matter, the symmetry energy per nucleon is a function 
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of the nucleon density. The form of this function has 
a bearing on the surface structure of finite nuclei, the 
structure of neutron stars, and the dynamics of heavy- 
ion reactions. This research was reviewed recently by Li, 
Chen, and Ko [6]. 

Because N — Z \s twice the eigenvalue Mt of the z 
component of the isospin T and most nuclear ground 
states are approximate eigenstates of with the eigen- 
value T{T +1) given by T = \Mt\, the symmetry energy 
may be conceived as the T-dependent part of the nu- 
clear ground state energy in the limit of isobaric invari- 
ance. Much theoretical and analytical work, reviewed in 
Sect. II, aims at describing its dependence on T. The 
empirical evidence seems, at least, compatible with the 
conjecture that, in a first approximation, the symmetry 
energy is proportional to the Casimir invariant of the 
isospin algebra, T{T + 1). 

In two previous brief articles [7], I discuss a schematic, 
micoscopic model which gives this T dependence approx- 
imately for low T. It is inspired by the Goswami's [8] 
observation that the potential of interaction, in a super- 
fluid nucleus, of the individual nucleons with the con- 
densate of Cooper [9] pairs, in the following referred to 
briefly as the "pair potential" , is not only nondiagonal 
in N and Z but also isobarically noninvariant. Frauen- 
dorf and Sheikh [10] point out that the symmetry energy 
may be conceived accordingly as an energy of rotation of 
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the condensate in isospace. An expression for the sym- 
metry energy proportional to T(T + 1) is then analogous 
to the well known expression for the energy levels of a 
quantal rotor. It should be stressed that it is, in this 
picture, the entire symmetry energy that is proportional 
to T(T + 1) and not only the contribution from the two- 
nucleon term in the Hamiltonian. In this respect the 
isorotational picture differs basically from several mod- 
els reviewed in Sect. II. For brevity, I call a rotation in 
isospace of a Cooper pair condensate "superfluid isorota- 
tion" . 

The model introduced in Rcf. [7] is designed to display 
this physics. It thus involves valence nucleons obeying 
a schematic Hamiltonian with a one- and a two-nucleon 
term which conserve the number A^, of valence nucleons 
and the isospin. To calculate states with arbitrary 
and T by a minimization, Lagrangian multiplier terms 
proportional to Av and Tz are subtracted from the Hamil- 
tonian, where A^ is the operator with the eigenvalues A^. 
This is equivalent to imposing neutron and proton chem- 
ical potentials. The Hamiltonian minus Lagrangian mul- 
tiplier terms is treated in the Hartree-Bogolyubov plus 
Random Phase approximation, which is known to lead 
to a separation of the collective degrees of freedom asso- 
ciated with spontaneously broken symmetries of a many- 
body system. Although the Hamiltonian is schematic, 
the principles of its treatment are thus general and may 
be applied to any energy functional of a Bogolyubov 
quasinucleon vacuum provided this hmctional is invariant 
under global gauge transformations and isobaric transfor- 
mations. 

The two-nucleon interaction of the model has a pairing 
and a particlc-holc part. When the latter is omitted, the 
theory is for a spherical nucleus equivalent to that of 
Ginoccio and Wesener [11]. Some of the results derived 
below are known from their study. However, the present 
methods arc quite different from theirs. The relation 
of the present work to that of Ginoccio and Wesener is 
discussed on the way in Sect. HI. 

In the present article, my model is analyzed to a much 
greater depth than in Ref. [7]. Its mathematical struc- 
ture is discussed in detail, its symmetries are explored, 
and the formulas used in the calculations are given ex- 
plicitly. Furthermore, the calculations have been con- 
siderably extended. In the idealized case of equidistant 
single- nucleon levels considered in Ref. [7] , the number of 
such levels has been enlarged by a factor more than 40 in 
order to rule out any spurious effect of the its finiteness. 
These calculations have been carried out for parameters 
appropriate for different mass numbers, and the results 
are described in a more general form than previously. A 
deviation from a linear T dependence of the contribution 
to the symmetry energy from the correlation energy cal- 
culated in the Random Phase Approximation (RPA) has 
been traced to the correlations of neutron-proton quasi- 
particle pairs and its origin understood. Finally, in order 
to approach a description of actual nuclei, calculations 
with Woods-Saxon single-nucleon levels have been car- 



ried out for several isobaric chains. 

To set up a background for the present work, I review 
in Sect. II previous theories of the symmetry energy dat- 
ing back to the era of the birth of nuclear physics in the 
1930ies. Certain aspects of some early models are elu- 
cidated in an appendix. The present empirical evidence 
as to the T dependence of the symmetry energy is also 
discussed in Sect. II. The formalism is then developed 
in Sect. III. and the calculations described and discussed 
in Sect. IV. After some brief remarks in Sect. V on the 
issue of an experimental signature of superfluid isorota- 
tion, the study is summarized in Sect. VI. 



II. THEORY AND PHENOMENOLOGY OF 
THE SYMMETRY ENERGY. A REVIEW 

Bethe and Bachcr [12] introduce the assiunption, which 
has since then been common in the literature, that the 
symmetry energy depends quadratically on N — Z. This 
means, in terms of isospin, that it is proportional to 
T^. A different T dependence is derived theoretically by 
Wigner [13]. He assumes that nucleons interact by a two- 
body force and makes the, now obsolete, assumption that 
this force is invariant under arbitrary unitary transforma- 
tions of the nucleonic spin and isospin. This implies, in 
particular, that any exchange force must be of the Ma- 
jorana type. Wigner then infers that the two-nucleon 
force gives a contribution to the symmetry energy equal 
to T{T + 4) times a factor which he supposes depends 
only weakly on T. If this factor is constant, the contri- 
bution of the two-nucleon force to the symmetry energy 
has a term linear in T besides the quadratic one. This 
is the celebrated "Wigner term". Wigner 's derivation is 
easily redone with the St/ (4) group of unitary transfor- 
mations of the nucleonic spin and isospin replaced with 
the SU{2) group of isobaric transformations. As shown 
in the appendix, a factor T{T+ 1) then replaces T{T+A). 
A term in the symmetry energy linear in T gives rise to 
a cusp at = Z in the curve of masses along an isobaric 
chain. Such cusps are found, in fact, in an analysis of 
measured masses by Myers and Swiatecki [14]. 

Wigner estimates the contribution of the nucleon ki- 
netic energy by the Fermi gas model, which gives a lead- 
ing term proportional to T^. In his model, the contribu- 
tions of the one- and two-nucleon terms in the Hamilto- 
nian thus depend differently on T. This is true for most 
models with a Hamiltonian composed of a one- and a 
two-nucleon term. 

In the framework of the spherical shell model, and as- 
suming conservation of isospin and seniority, Talmi and 
Unna [15] show that, for nuck;i whose valence rniclcons 
occupy a single j-shcll, the symmetry energy is propor- 
tional to T{T+1). In the spherical shell model, the single- 
nucleon term in the Hamiltonian has the role played in 
the Fermi gas model by the nucleon kinetic energy. When 
all valence nucleons occupy the same j-shell, the sum 
of their single-nucleon energies is constant for a given 
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mass number, so the entire symmetry energy stems from 
the residual two-nucleon interaction. While the result of 
Talmi and Unna thus agrees superficially with that of 
Wigner's argument applied to the isobaric SU{2), it is 
shown in the appendix to imply that the basic assump- 
tion of the latter, namely that the average interaction 
energy in two-nucleon states with a definite symmetry in 
position and spin is independent of T, does not hold for 
the isobarically invariant isovector pairing force acting in 
a single j-shell. 

In isobarically invariant shell model calculations for 
A = 24 and A = 48, Satula et al. [16] find that omitting 
the interaction of isoscalar nucleon pairs essentially elimi- 
nates the deviation from a quadratic T dependence of the 
symmetry energy. Since the ^ = 48 nuclei whose calcu- 
lated binding energies are analysed in Ref. [16] belong 
to the 1/7/2 shell, this seems to indicate, in view of the 
result of Talmi and Unna, that the interaction of isovec- 
tor pairs used in the calculation has a major component 
which does not conserve seniority. In line with an earlier 
study by Brenner et al. [17], Satula et al. also consider a 
certain linear combination of the binding energies of nu- 
clei differing in TV and Z by small numbers. These linear 
combination are constructed so as to filter out a posi- 
ble enhancement of the binding energy for N = Z. Both 
Brenner et al. and Satula et al. find that, for a range 
of sd and fp shell nuclei, their linear combinations van- 
ish essentially when the interaction of isoscalar pairs is 
turned off in their shell model calculations. They hence 
infer that the observed deviation from a quadratic T de- 
pendence of the symmetry energy stems from this part 
of the two-nucleon interaction. 

Myers and Swiatecki [14] fit the measured nuclear 
masses with a formula where the deviation from a 
quadratic dependence of the symmetry energy on N — Z 
decreases exponentially with \N ~ Z\. More recently, 
Myers [18] suggests that a term in the symmetry en- 
ergy linear in | A'' — Z\ could arise because two nucleons in 
identical states of orbital motion are more strongly bound 
than other nucleon pairs. He assumes that the single- 
nucleon levels are fourfold degenerate and the states of 
each level differ only by the directions of their spins 
and isospins. Thus he neglects the single-nuclcon spin- 
orbit potential. His independent-nucleon model has, in 
fact, Wigner's 5f (4) symmetry. When N neutrons and 
Z protons occupy the lowest levels in a potential well, 
and N and Z arc even, the number of pairs in identical 
states of orbital motion is then equal to 3A/2 —\N — Z\. 
In the presense of a single-nucleon spin-orbit potential, 
this is the number of pairs of nucleons in identical or 
time reversed states of orbital motion and spin. Jensen, 
Hansen, and Jonson [19] work out a version of Myers's 
argument which maintains the S'C/(4) symmetry of the 
independent-nucleon model but includes exchange terms 
in the assumed delta-force interaction of the nucleons. 
In Ref. [20], Myers and Swiatecki explain the extra bind- 
ing of nucleon pairs in identical states of orbital motion 
by the congruence of the nodal surfaces of the single- 



nucleon wave functions. They call the resulting term in 
the nuclear binding energy, accordingly, a "congruence" 
energy. In Ref. [21], these authors consider only pairs of 
a neutron and a proton and give for the number of par- 
ticularly strongly bound pairs of this kind the expression 
{A—\N — Z\)/2, which is the number of neutron-proton 
pairs in identical states of orbital motion and spin. Since 
a congruence of nodal surfaces requires that the station- 
ary states of a nucleon have unique orbital wave func- 
tions, and thus that the single-nucleon spin-orbit poten- 
tial is neglected, counting only nucleon pairs with parallel 
spins seems, however, inconsistent with the congruence 
picture. 

A contribution to the symmetry energy proportional 

to T(r+ 1) is pointed out by Bohr and Mottelson [3] to 
arise from the separable particle-hole interaction which, 
in the Hartree approximation, generates a term in the 
single-nucleon potential proportional to T • t, where t is 
the single-nucleon isospin. Satula and Wyss [22] study 
the symmetry energy in Hartree-Fock calculations with 
Skyrme forces where they ensure isobaric invariance by 
omitting the Coulomb force and assuming equal neutron 
and proton masses. They find that when the isospin- 
dependent parts of the Skyrme forces are omitted, the 
symmetry energy is roughly proportional to with T 
taken here equal to Mt. The isospin-dependent parts 
give additional contributions which are nearly propor- 
tional to T(T-f- 1). They thus act similarly to the sep- 
arable particle-hole interaction of Bohr and Mottelson. 
The quadratic T dependence of the remainder of the sym- 
metry energy may be understood to result from a redis- 
tribution of the nucleons on their self-consistent energy 
levels. Fig. 4 of Ref. [22] shows that when a Bardeen- 
Cooper-Schrieffer [23] pairing term is added to the en- 
ergy functional, the contribution of the isospin-dependent 
forces to the symmetry energy is no longer proportional 
to T(T + 1), where T is now taken equal to (Tj). A 
roughly quadratic dependence on T can be inferred from 
the article's Figs. 2 and 4. This is explained in Sect. IV B 
below. 

In a model with nucleons in a deformed potential and 
a pairing force acting on both isoscalar and isovector nu- 
cleon pairs, Satula and Wyss [24] find that an approx- 
imate particle number projection allows isoscalar and 
isovector components of the pair potential to coexist. In 
calculations for doubly even nuclei in the Jp shell, the 
isoscalar component turns out to be for certain values 
of the coupling constants particularly large for N x Z. 
which leads to an enhanced binding of such nuclei. The 
authors therefore suggest that the Wigner term arises 
from isoscalar pairing. In Ref. [25], they consider an 
isoscalar pair potential with a different structure. With- 
out particle number projection, Civitarese, Reboiro, and 
Vogel [26] obtain a similar enhanced binding for N k Z 
with an isovector pairing force which acts with a larger 
strength on neutron-proton than on neutron and proton 
pairs and is thus isobarically noninvariant. 

The empirical evidence as to the T dependence of the 
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symmetry energy is ambiguous. Thus in the recent analy- 
ses by Royer and Gautier [27], Royer [28], Kirson [29, 30], 
and Dieperink and Van Isacker [31], several terms which 
are plausible parts of a mass formula compete to improve 
the fit to the empirical masses. It seems safe, though, to 
conclude that a Wigner term in some form is called for 
and that a symmetry energy propotional to T(T + 1) is 
compatible with the data. If a factor T{T+x) is assumed, 
the tendency is that x is somewhat less than one. The 
analyses by Zeldes [32] and Kirson [30] also indicate that 
the deviation from a quadratic T dependence is not con- 
fined to the region of nuclei with approximately equal A'^ 
and Z, such as implied by the exponential parametriza- 
tion of Myers and Swiatecki [14] and the theories of Bren- 
ner et al. [17], Satula et al. [16], Satula and Wyss [24, 25], 
and Civitarese, Reboiro, and Vogel [26]. 



III. THEORY 

A. Hamiltonian 

The Hamiltonian considered may be thought of as that 
of a spherical or deformed shell model. It is written 

^ = Y^ ^A^3 + 1 ^Jfcima]4«m«i > (1) 

j jklm 

where aj annihilates a nucleon in the state |j) and 

the summations runs over an orthonormal set of single- 
nucleon states spanning a valence space. The basic 
single- nucleon states form quadruples \j) = Iqctt) with 
a common energy ej = eg. where qa is q or ^, and r 
is n, denoting a neutron state, or p, denoting a pro- 
ton state. These states are related by j^r) = \qT) and 
jgp) = f-jgn), where a bar over a ket denotes time re- 
versal, and t± = ± ity in terms of the single-nucleon 
isospin t = {txTty,tz). Note that time reversal is defined 
so as to commute with t-, so tx and tz are even and ty 
odd under time reversal. 

The interaction matrix element Vjkim has a pairing and 
a particle-hole part, which are treated differently below, 

„ _ ,,Pair , ph 

The pairing part is the matrix element of the isobarically 
invariant isovector pairing force, 

v^^l = -G{jm.{m\t\l}, \j) = -2ity-\j). (2) 
The components of 

P = {Px, Py, P.) = 71 E iMMJKaj 

jk 



are given by 

P+ = Px + iPy = «3pV = ' 

? 

P-=Px- iPy = -V2Y,a^n%n = '^Pn , (3) 
9 

Q 

It follows from = —{k\t\j), which is easily verified, 

that P commutes as an isovector with the isospin 

T = {Tx,Ty,Tz) =J^{j\t\k)a]a^. 

jk 

Note 

T^=Ti=Tx + iTy = at^„a,,p , 

q(7 ^ ^ 

q<T 

The particle-hole part of Vjkim is the matrix element of 
the separable interaction of Bohr and Mottelson [3] , given 
by 

v%m = '^{jm-{mm) (5) 

with a coupling constant k. I call this interaction the 
"symmetry force" . 
The Hamiltonian (1) commutes with T and 

i 

It therefore has a complete orthogonal set of eigenstates 
which are also eigenstates of Ay, T^, and T^. I denote 
the eigenvalues of these operators hy Ay, T(T + 1), and 
Mt- The numbers Nr of valence neutrons and protons 
arc then given by Nr = + 2miMT with in, = ±1/2 

for r = n and p. Since we are concerned with isobaric 
multiplets where these numbers are even for Mt = T, 
it follows that is even, T is an integer, and Av/2 
and T have equal parities. The eigenstates of H form 
degenerate multiplets with Mt = T,T — 1, . . . , —T. 

If G = 0, the lowest eigenvalue of H for given A^ and 
T is 

^ = 2^ E e, + ^n(T{T+l)-lAy), (6) 

where the Fermi level is such that Nr/2 levels sat- 
isfy Cq < Xr with NT = Ay/2 + 2mtT. The eigenvalue 
of the symmetry force in Eq. (6) reflects X^a^b ta ■ tb = 
— *a =3/4, where ta is the isospin of the 

ath nucleon. The energy E given by Eq. (6) is a convex 
function of A^ and T and an increasing function of T for 
given Ay. These properties of the lowest eigenvalue of H 
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are likely to persist for G 7^ 0. Since the actual ground 
state energy, reduced for the electrostatic energy, of a 
doubly even nucleus with N = Z is not a convex function 
of A, the Hamiltonian (1) therefore does not reproduce 
the absolute reduced energies, but it may reproduce their 
differences for a fixed A. 

This Hamiltonian is evidently very schematic. Its rel- 
evance for a study of superfluid isorotation is due to the 
fact that it obeys the symmetries that are violated by 
the pair potential of Cooper pairing, global gauge in- 
variance and isobaric invariance. Another virtue of the 
Hamiltonian is its simplicity, which makes its behavior 
transparent. Other symmetries of a more realistic nu- 
clear Hamiltonian than the global gauge invariance and 
the isobaric invariance are violated by the present one. 
Like every shell model Hamiltonian, it is thus not invari- 
ant under translational and Galilean transformations. If 
the single-nucleon energies are derived from a deformed 
potential, as in the calculations in Sect. IV C below, it 
is also not invariant under rotations in space. If the 
Hamiltonian would obey these symmetries, one could 
have solutions of the Hartree-Bogolyubov problem which 
would break the symmetries and thus give rise to Nambu- 
Goldstone solutions of the RPA problem additional to 
those discussed in Sect. IIIH below. The translational 
and Galilean invariances are necessarily broken by the 
solution of the Hartree-Bogolyubov problem. Imposing, 
in particular, rotational invariance would be important 
in a study of rotationally excited states. For the present 
study, which deals with nuclear ground states, and where 
excited states arc considered only in so far as they are iso- 
baric analogs of ground states, rotational invariance may 
be assumed to be less important. 



angular momentum, a "Routhian" , a term borrowed from 
analytical mechanics. 

C. QuEisinucleon vacuum 

I now set out to calculate R approximately by per- 
turbation theory starting from a vacuum |$) of Bo- 
golyubov [34] quasinucleons. The state |$) is determined 
up to a phase by 

a.|$)=0 (7) 
in terms of a complete set of quasinucleon annihilators 

k 

obeying 

{a^.,aj=0, {aj,al}=djk- 

It will be assumed to minimize the "Hartree-Bogolyubov" 
Routhian Rub given by 

Rhb = Ehb - A(iv) - , (8) 

^HB=$^e,(a]a.) 

3 

jklm 

= {Ho)-G\{P)f + ^K{T)\ (9) 



Routhian 



The eigenstates of H with the quantum numbers Ay, 
T, and Mt are also eigenstates of 



R = H-\Ay 



where A and fi are parameters, and the lowest eigenvalues 
E and RoiH and R for given A^, T, and Mt are related 
by 

R = E- XAy- iiMt . 

If _E is a convex function of A^ and T and an increasing 
function of T for given Ay, then for any set of A^ and T 
there exists a set of A and /z with /li = for T = and 
> for T > so that R is minimal for this set of A^ 
and T and Mt = T. Hence the lowest eigenstates of H 
for given A^ and T and with Mt = T is also the lowest 
eigcnstate of R for some values of A and /i > 0. Because 
the allowed values of Ay and T form a discrete set, A 
and are not unique functions of these variables. They 
can be chosen freely within certain limits. Pashkevich 
and I [33] call a quantity analogous to R, involving the 



In these expressions, the expectation values are in the 
state 1$). This is my convention from now on unless 
otherwise stated in the context. If Rub is minimal, |$) 
is an eigenstate of 



X - dRuB 
Ho - G{{P)* ■ P + (P) ■ P^) 

- XAy - [iTz , 



k(T) • T 



(10) 



where the convention has been introduced that a Greek 
letter subscript takes a value j or with a-^i = '^]- 
(A similar notation is used in an early article [35] by 
Vogel and me.) The partial derivatives in Eq. (10) refer 
to the expressions (8) and (9) with all (a^a^) considered 
mutually independent. Since |$) is an eigenstate of i?o, 
the annihilators aj can be chosen as eigenvectors of the 
linear map a 1— > [a, i?o]. 

If the isovector (T) is different from the zero vector, it 
can be rotated into an arbitrary direction by an isobaric 
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transformation of | $) . This does not change the first two 
terms in the expression (8). For /x > 0, the isovcctor (T) 
therefore points into the z direction if i?HB is minimal. 
For = 0, the value of JRhb does not depend on the direc- 
tion of (T), which can therefore be just assumed to point 
into the z direction. Generally, then, (Tx) = (Ty) = 0, 
which also holds if (T) is the zero vector. 
^ The exact ground state of R is an eigenstate of 
Nr = Ay/2 + 2mtTz with even eigenvalues and therefore 
an eigenstate of e~'^'^^^ with the eigenvalues one. If 1$) is 
an eigenstate of e~^'^^^ , the eigenvalues should therefore 
be one as well; that is. 



(11) 



Suppose i?HB has been minimized with these con- 
straints, and consider an infinitesimal variation \5^) 
which violates the constraints. Any quasinuclcon 
vacuum is an eigenstate of e"""^^ = e^*'^^"e^™^p, 
which has the eigenvalues ±1. By continuity, there- 
fore, e-"^°e-*'^*p|(5$) = Since e"*'^*- has 
the eigenvalues ±1, the variation must then satisfy 
g-iwN^ |i5<[)) = -|(5$). It is easily verified that Rub is sta- 
tionary with respect to such a variation, and it will be 
seen in Sect. IIIH that i?HB is locally minimized by a 
state which satifies the constraints (11). I cannot prove 
that it is also globally minimized by this state; this is- 
sue is discussed a little further in Sect. IIIH. Anyway, 
I impose from now on the constraints (11). This entails 
(Pz) = 0, so Eq. (10) becomes 

Ro = Ho-J2{A*P,+A,P^) 

T 

+ (k(T,) - ^) • - Aiv , A, = G(P,). 

A transformation |$) i— » e~^^-r'^^^^\^) with suitable 
angles 4>r makes > 0. This leads to an expression for 
i?o in the form of the single-quasinucleon Hamiltonian 
of the theory of Cooper pairing of Bogolyubov [36] and 
Valatin [37], 

Ro=Ho- ^(A,(P, + Pi) + XrNr) , 

r 

X^ = X + mt{n-K{T,)), (12) 
so that the solution of the eigenproblem 

[aj,Ro]= E^aj , > , 
is obtained immediately from this theory. 




E„. 



= E-,r = ^{e,-KY + A%. (15) 



Because is an eigenvector of a i-^- [a, i?o] with the 

eigenvalue —E-, either or a] could annihilate |$). It 
is shown, however, in Sect. IIIH that a quasinucleon vac- 
uum annihilated by an operator cannot locally min- 
imize i?HB if it satisfies Eq. (11). The operators (13) 
therefore annihilate |$), and Eq. (11) holds by the con- 
struction of 1$). 

The relations A^ = G{Pr) are equivalent to 



0. 



(16) 



Unless otherwise stated in the context. A,- > will be as- 
sumed. It remains to fix A,-. The conventional way of the 
theory of Cooper pairing, which I shall follow, consists in 
demanding 



(TV.) = 2 ^ vl, = Nr = ^ + 2m,T . (17) 
g 

This implies 



E = EuB + {R— Rhb) ■ 

In other words, corrections to Ehb can be calculated as 
corrections to Rub- From Eqs. (12) and (17) one gets 

At- = A + mj(/i — kT) , 
A=(A„ + Ap)/2, n = Xn-Xp + KT, (18) 



and Eq. (9) becomes 

EuB = Eo + £^pair + ^I^T"^ , 

Eo = {H^) = 2Y,vl,e^, 



(19) 
(20) 
(21) 

Eqp. Furthermore, An — Ap and, hence, /U 



E, 



2 A2 + A2 

pa. = -G\{P}f = 



If and only if T = 0, one has Uqn = u 
and Eqj 

are increasing functions of T, so T = 0, An = Ap, and 
/i = arc equivalent. For T = 0, also An = Ap, so in 
terms of its Cartesian components, (P) is purely imagi- 
nary and points into the y direction. Other states with 
the same i?HB are generated in this case by arbitrary 
global gauge transformations and isobaric transforma- 
tions of 1$). These transformations generally change the 
complex argument and the direction of (P) . Rotating, in 
particular, (P) into the z direction by an isobaric trans- 
formation yields An = Ap = and Anp = G{Pz) ^ 0. 
This solution of the Hartree-Bogolyubov problem for the 
isobarically invariant isovector pairing force is mentioned 
by Engcl et al. [38] as alternative to the former. Both 
solutions are in fact just two out of an infinity of equiv- 
alent solutions related by global gauge transformations 
and isobaric transformations. In general, these solutions 
violate the isobarically noninvariant constraints (11). 
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D. RPA 



by 



In my brief articles [7], I derive the RPA part of the 
theory from a boson expansion. This allows some short- 
cuts by reference to the literature. Boson expansions are, 
however, ambiguous due to the noncommutability of bo- 
son field operators. Following Thoulcss [39], I base the 
following derivation on perturbation theory. 

For any Hamiltonian Ti, and any pair of operators X 
and Y, I define the propagator 



gt{X,Y,t,'H) = -i{T{X'{t)Y'm) , 
X{t) = e'^*Xe-*^*, X' = X-{X), 



(22) 



where T{...} indicates time ordering and the expec- 
tation values are in the ground state of Ti. It is eas- 
ily proved by Thouless's [40] method that gt{X,Y,t,n) 
is the sum of all Fcynman diagrams, referring to some 
independent-particle Hamiltonian, where the vertices 
representing X are linked to the vertices representing Y 
and with no unlinked part. In particular, the subtraction 
of {X) in the last of the equations (22) cancels diagrams 
without an unlinked part where the vertices representing 
X are not linked to the vertices representing Y. li X 
and Y are linear combinations of products of even num- 
bers of fcrmion field operators, the Fourier transform of 

gt{x,Y,t,n) is 

/oo 
e"^%{X,Y,t,n)dt 
-OO 



= X' 



(J — {H — £ — irj) 



Y' -Y' 



• + {n-e- in) 



X' 



where £ is the lowest eigenvalue of H and > is in- 
finitesimal. I consider propagators 

G{X,Y,u)^g^{X,Y,uj,R), 
GoiX,Y,uj) = gUX,Y,LO,Ro), 

and Fcynman diagrams referring to Rq as the 
independent-particle Hamiltonian. In particular, 



Go{ai^ai^,a^a^,u}) = otherwise. 

The propagator G{a^ai^,a^a^,(jj) is approximated by 
the sum of the diagrams 



(23) 



where a "lens" represents Ga{ai^a^,as^a^,uj), and a 
dashed line the interaction matrix element w^^ given 



''^j-H,k--'m ~ '"jklm ' 

"'t«,AM = Otherwise. 



(24) 



Due to the separable form of these matrix elements ac- 
cording to Eqs. (2) and (5), the diagrams (23) involve 
free propagators Go {X, Y, cu) with many coherent terms 
over the configurations of two quasinucleons. They there- 
fore give a large contribution to the exact propagator 
G{a^a^,a^^aij^,u>). Summation of these diagrams results 
in the RPA equation 

G{a^a^,ay^a uj) = Go{a^a^,a^a u) 



TTpar 



In a basis of a complete set of independent operators 

af,aj and their Hermitian conjugates a^._iaj^_i, this 
equation takes the matrix form 

G(a;) = Go(w) + Go(a;)VG(w) = (Go(w)-^ - V)-i , (25) 

where 

V.«,Am = Yl K«p).K(«-«r)AM^-P,-r (26) 
npar 

with 

iX).. = {{^W.,X}). (27) 
Eq. (25) implies that G{lo) has poles where 
Go(a;)-^-V 

is singular. 

It is convenient to introduce also block matrices corre- 
sponding to a division of the basis into a first part con- 
sisting of the operators ctf-aj and a second part consiting 
of the operators a ._ia,_i with a common order of the 
pairs (j, k) in both part. Since 

(Go(a;)-i) = a; - {E. + E, - irj) , 

(Go(w)-')^^^^ = otherwise, 

the poles of G(w) arc then the eigenvalues of 

R=(_°^ ;)(E + V), (28) 

where and 1 denote the zero and identity matrices and 

^j--k--,kj = = E^ +E,-iv, (29) 

E.^,A^ = otherwise. 
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The properties of R and 

K=(; J)(E + V)=(J _°,)R (30) 

are discussed by Thouless [39]. I shall follow his discus- 
sion partially. Because E + V is symmetric, 

R" = (E+v)(; V)=-(-°i V)- 

Therefore, R has an eigenvalue —w for every eigenvalue 
u. For 77 = 0, K is Hermitian. If x is a right eigenvector 
of R with the eigenvalue w, therefore 

xt(J _ojR = xtK 
= xtKt=xtRt(J _°,)=.*xt(j _o^), (31) 

so x^ is a left eigenvector with the eigenvalue uj*. 

Two right eigenvectors x^ and X2 with the eigenvalues 0;^ 
and a»2 satisfy 

so 

xl(J _°i)x2 = if ojI^uj^. (32) 
If K is also positive definite, the matrix 

kV2rk-V2 = kV2(J _o^)kV2 

is Hermitian, so the eigenvalues of R are real, and its 
right eigenvectors span the space of cohimn matrices of 
their dimension. (This simple argument is due to Ring 
and Schuck [41]; Thouless [39] has a more complicated 
one.) If X is a right eigenvector of R with the eigenvalue 
UJ, one then gets from Eq. (31) 

wxt^J -i)x =x^Kx >0, (33) 

so Lo cannot be zero. When the imaginary term in 
Eq. (29) is included, R acquires the additional term 

5R = -i,(j 

It follows from Eq. (31) that the resulting change 5u of 
ijj is given by 

5a;xt(j _o^)x =xt(j _o^)jRx =-^r;|x|^ 

A comparison with Eq. (33) then shows that the positive 
eigenvalues of R move into the lower imaginary halfplane. 

It will be seen in Sect. IIIG that K is, actually, only 
positive semidefinite for 77 = 0. To escape the complica- 
tions arising from the eigenvalues zero of K in this limit, 
it is understood in the next section that an infinitesimal 
positive definite term has been added to K. 



E. Correction R — Run 

The correction R — i?HB can be divided into two parts, 

R-Rhb = {{R) - Rub) + {R- (R)) ■ 

A general result derived by Goldstone [42] implies that 
the second bracket in this expression is tlic^ sum of all 
linked Feynman diagrams without external lines and with 
at least two interaction lines. The first bracket is the part 
of the expectation value of the two-nucleon interaction in 
the quasinucleon vacuum that is not included in i?HB • It 
is given by 

(R) — Rhb 

jklm 

= i 5Z \^^Jfcim({«]4 - {a]al),a^ai - (a„a;)} 

jklm 

+ [a]4 

+ at[4a,„,a;])y 
It will be seen that the part 

jklm 

+ (4"( ~ ("-hi )) (4am - (4«m)) ) (34) 

of these terms can be combined with terms in i? — (^?). 
The remainder c can be calculated from Eqs. (2), (5), 
and (17), 

^=1 Z!(^^jfc;m[44'amai] + ^^jfeJm4K«m>ai]) 
jklm 

= f(G(2rf-^v)-|KAv), (35) 

where M is the dimension of the valence space. The 
factor 3/4 is just t^. It is seen that c docs not depend 
on T and thus does not contribute to the symmetry en- 
ergy. The pairing-force part of c vanishes when the va- 
lence space is halfway filled. The symmetry-force part is 
5K;( J2a^b ^a'ib — T^) i Compare the remark after eq. (6). 

The expression (34) is equal to the expression (37) be- 
low with n = 1, which is the first diagram in the series 
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I approximate R — i?HB ^ c by the sum -Erpa of this se- 
ries. This is expected to be a good approximation for the 
reason mentioned in connection with the diagrams (23). 
The nth diagram equals 

where the denominator 2n appears because the diagram 
has 2n equivalent vertices. Hence 



-E^RPA = 



.f:itr(VGoH)" 



n=l 



dw 
2^ 



(38) 



When the positive eigenvalues of R for 77 = are denoted 
by LOn, and the set of pairs (j, k) of subscripts of the basic 
operators ataj by S, the integrand in Eq. (38) can be 
expressed as follows. 

-$:itr(VGo(o.))" 

n=l 

= trlog(l - VGo(w)) = logdet(l - VGo(a;)) 
= logdet((Go(u;)-^-V)Go(a;)) 

= logdet((a;-R)(a;-(_°^ J)e)"'^ 



log! J|(w-(w„-ir?))C 

\ n 



I U,k)es J 
= - (w„ - ir])) + log{u> + (a;„ - irj))^ 

n 

(i,fe)65 

+ \og{u + {E^+E^^-iij))^ . 

Since this is proportional to uj~'^ for large cj, the inte- 
gration path can be extended with an infinite semicircle 
in the lower imaginary halfplane. When this semicircle 
is deformed into a linear path running backward below 
the real axis at the distance 2ry, the integral collects con- 
tributions only from the discontinuity of log z for z < 
in the terms in the integrand which have this cut in the 
lower imaginary halfplane. The result is 



-E'rpa = \ E'^" 



iE,+E, 



(39) 



This expression has a very simple interpretation if lo^ 
and Ej + Ej. are conceived as frequencies of harmonic 
oscillators. Then, Srpa is just the change in total oscil- 
lator zero point energy induced by the interaction of the 
quasinucleons. 



Truncating the expansion in Feynman diagrams to the 
sum of the diagrams where a pair of quasinucleons cre- 
ated together is also annihilated together is ccjuivalcnt 
to treating such a pair as a single boson. Therefore, 
the result (39) is the same as obtained in the quasi- 
boson approximation. To my knowledge, the deriva- 
tion above does not appear in the literature. A some- 
what different derivation by contour integration is due 
to Shimizu et al. [43]. In the quasiboson approximation, 
Marshalek [44] derives an expression for the RPA correc- 
tion to a Hartree-Bogolyubov energy which corresponds 
to the present c -I- -Brpa when the difference between the 
Routhians considered is taken into account. 



F. Matrices E and V 

Prom now on, ry is set to zero without further notice, 
so E is real. For any operator X, I denote by the corre- 
sponding sans-serif symbol X the column matrix of the 
brackets (X) defined by Eq. (27). From Eqs. (2), (5), 
(24), and (26fthen follows 



-G 



((;; 



p* + p p 



X? 0)) 



+«(i(T+T^+T_T^)+T,T^ 



(40) 



It is seen from Eqs. (3), (4), (13), (14), and (27) that 

P^, P^, T_|_, and T^ are real. Furthermore, if X is any 
of these matrices, then, with j = qar and k = q'a'r', 
its elements and (X)^._i^_i are nonzero only for 

q' ^ q and a' 7^ a. Therefore, only within the subspace 
spanned by the corresponding basic operators do the 
eigenvalues of R differ from those of (^^^ J^E and thus 

contribute to the difference in Eq. (39). Within this sub- 
space, one can choose j = qr and k = qr'. 1 call the space 
spanned by such basic operators with t = t' the rr space, 
that spanned by such basic operators with t ^ t' the np 
space, and the direct product of the nn and pp spaces the 
nn+pp space. Since for any of the aforementioned ma- 
trices X, the corresponding operator X is time reversal 
even, it follows from Eq. (27) that {X)_ = (X)_ 

' V / \ / qpqn \ / qnqp 

and (X)(^^)_,(_p^_, = (X)(^p)_,(_^^_,. Therefore, if the 
basic operators of the np space are replaced with their 
linear combinations 



qq± V2 



72Hn4p±4p4n) 



[atat] 



(41) 
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only the corresponding brackets 



G. Deflniteness of K 



(X)_.=^((X)_ ±(x)_ 

V y qq± v2 \^ / qpqn V / qrii 
Wg-ig-ii = ;^((^)(qn)-i(gp)-i ^ ('^) (gp)-i(9n)-0 

with the subscript + are different from zero. The only 
nonzero elements of E in the basis of the operators (41) 

are 



Eg-ig-i± ,qq± — Eqq± ,g-iq-i± — -Bgn + -E'qp • (42) 



For the calculation of E^rpa one therefore needs to keep 
only the part of the np space spanned by the operators 
with the subscript +. From now on, I call this subspace 
the np space. 

Eqs. (3), (4), (13), and (27) give 



qq+ / "ijn^qp ; 



(-^-r)(9T)-MqT)-i ~ y^qT^^qrPT 

{P.)q-.q-.^ = {[^^]qq^P. 



q-^q 



I a''' a'''] 



qq+ 



Cr-)-qq. = {T.)q-^-q-^. = {T-[-'-' 



' q-^q 



(43) 



i^^)qTqT = ("^^)(<}r)-l(9T)-l = 



vt 



[X)^ = otherwise, 

where a is qrqr, {qT)~^ {qT)~^ , qq+, or q~^q'^^+. The 
elements of V connecting the nn+pp and np spaces are 
seen to vanish. The part of V acting in the nn+pp space 
is given by the terms in the sum (40) with P^, P^., and 
T^, and the part acting in the np space by the terms with 
P^, Pt, andT±. The nonzero elements of E in the nn+pp 
space are 



.grgr " ^grgr , (qT)-i(gT)-i " 2^^^ . (44) 



In the np space, they are given by Eq. (42). Note that 
for K = 0, the elements of V connecting the nn and pp 
spaces vanish. When also 1-1 = 0, the nn, pp, and np parts 
of E or V are identical. 



Since the matrices E and V are real, the symmetries 
mentioned in Sect. HID imply 



E + V=(B A), 

where A and B are symmetric. Hence follow 
K = R 

whence 



(45) 



R' 



A+ B 

A - 

A - 

A + B 



(46) 



I shall prove that K is positive semidefinite. 
Prom T = rt and T, = tI follow T=(° i'lT, and 



1 

_|_ - luiiuw ' z - oj 

I oj'^^- Since k is positive, the symmetry force 
is seen from Eqs. (30) and (40) to contribute, then, a 
positive semidefinite term to K. It is therefore sufficient 
to consider the case k = 0. From the remark at the end 
of Sect. IIIF, it follows that, in this case, the nn and 
pp spaces can be treated separately. Ginoccio and We- 
sener [11] prove that K is positive semidefinite in these 
spaces for k = 0. The present proof is somewhat different 
from theirs and extended to the np space. It also gives 
the dimension of K's kernel. 

With the tt' parts of A and B denoted by a-rr' and 
b^^/, it is according to Eq. (46) sufficients to prove that 
the matrices Stt'^^tt' are positive semidefinite for k = 0. 
Eqs. (40), (42), (43), (44), and (45) give 



Gp 



rr'=p"TT'=p 



(47) 



where Ct-t' is the diagonal matrix with the diagonal ele- 
ments + Eqr' and Ptt'± the column matrix with the 
elements 

PqTT'± = UqrUqr' ± VqrVqr' . (48) 

For any column matrix x, thus 

= (1 - G/,,.^)xte,,,x 

with 



/• 



tt'± 



Prr'±« 



'tt' "rr'ih 



E 



PqTT'± 



Ear + E, 



qr' 



where equality holds if and only if x is proportional to 
&-l.Prr+- Eqs. (16) and (48) give 



rr'ib 



< 



E 



Ea 



Ea 



E, 



+ 



•qr 



E. 



qr' 



1 



11 



where, as it follows from the remark at the end of 
Sect. IIIC, both equalities hold if and only if the sign 
is + and either t — t' or fi — 0. Hence, the matri- 
ces Btt' ± bxT' are positive semidefinite. Only a matrix 
aT-T-/ — brr' can have an eigenvalue zero, which it has if 
either t = t' or /i = 0. The multipUcity of these eigen- 
values is one. They are discussed from another point of 
view in Sect. IIIH. 

Let the symmetry force be included again. Since, as it 
has now been shown, A + B is positive definite, one can 
apply to R one more similarity transformation. 



of the positive eigenvalues of M. The eigenvalues of M 
are calculated most efficiently from 



M' = (A + By/^M{A + B)-i/2 = (A + B)(A - B) . (49) 



That the dimension of the matrix to be diagonalized can 

be reduced in this way when R is real is pointed out 
by Ring and Schuck [41] in the case when K is positive 
definite. 



R" 



(A + B)^ 




(A + B) 



-1/2 



R' (A 



(A + B) 



1/2 



(A-B)(A + B) 




1/2 



1/2 

(A+B) 

) = (?X)- 



1/2 



The matrix M is real, symmetric, and positive semidef- 
inite. It therefore has only nonncgativc eigenvalues tj^, 
and the corresponding eigenvectors x span the space of 
column matrices of their dimension. If > 0, then 



(UJX 



are eigenvectors of R" with the eigenvalues uj. Note that 
LO can have both signs, so two eigenvalues of R" with 
opposite signs correspond to each such eigenvalue of M. 



If a; = 0, one has 



H. Nambu-Goldstone solutions 



It follows from a theorem proved by Thouless [45] that 
a general Bogolyubov quasinucleon vacuum in the vicin- 
ity of |<E>) is obtained from |$) by a transformation of the 
form 



|$)^^exp| Cfc,"Hll^)' 
k(j,fe)es 



(50) 



R' 



fffO 



0, R'' 



The column matrix in the first of these equations is an 
eigenvector of R" with the eigenvalue zero, whereas the 
first one in the second equation is not an eigenvector of 
R". Following Thouless [39, 45] and Marshalek and We- 
scner [46], one can interpret the linear combinations of 
two-quasinucleon annihilators and creators correspond- 
ing to these column matrices as proportional to leading 
terms in expansions of a conserved momentum and its 
conjugate coordinate. Together, all the pairs of column 
matrices mentioned above span the space of column ma- 
trices of their dimension. 

It follows from Eq. (30) that K and R have a com- 
mon right kernel, whose dimension I denote by dk. For 
K = 0, the preceding discussion implies dk = 2 if > 
and rfk = 3 if /i = 0. Since rfk is also the dimension of the 
right kernel of M and 

det(a; - R) = det(a; - R") = det(a;^ - M) , 

the characteristic root zero of R has the multiplicity 2(ik. 
When K is made positive definite by an infinitesimal per- 
turbation, as it was assumed to have been done for the 
purpose of the derivation in Sect. HIE, the eigenvalue 
zero of R then splits into dk pairs of nonzero infinitesimal 
real eigenvalues with opposite signs. Being infinitesimal, 
these eigenvalues do not contribute to E'rpa, which can 
therefore be calculated from the positive eigenvalues of 
the original singular matrix R. These are the square roots 



where C^,^ are complex parameters, and A/" is a normaliza- 
tion factor. Since expectation values in the transformed 
vacuum do not depend on its phase, they are functions 
of (^f.j . Thouless's results in Ref. [45] imply that E + V 
is the Hessian matrix of i?HB with respect to (j^j and 
~ Cki ^'^^ Cki ~ ^ and that, if some symmetry of 



EuB is violated by the quasinucleon vacuum, special so- 
lutions of the eigenproblem of R result from considering 
variations of ]$) within the symmetry group. These solu- 
tions are analogous to the Nambu-Goldstone [47] boson 
solutions in field theories with a so-called "spontaneous 
symmetry breaking", that is, a violation of a symmetry 
by the vacuum. 

In the present case, the variations to be considered are 



\S^) = -iSxX\^)., 



(51) 



where X is any one of the Hermitian operators N^, 
and Ty, and 5x is infinitesimal and real. Because the first 
two terms in the expression (8) are invariant under global 
gauge transformations and isobaric transformations. 



5R 



HB 



5X 



■ ili{Y), Y=[T,,X]. 
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whence for <^i^j = 0, by S-Rhb/^Cik — follows 



E 



(A-i,/*-!) e 5 



E 



(fi. A) G 5 or 



d [ORbbK 



E 

{fi, X) € S or 

(A-\/x-i) e s 

d 6RnB . d{Y} . 

= (52) 



5C« -^x ■ 5C, 

For Cfc^. = 0, Eqs. (50) and (51) give 

5C._,,_, = SCkj =i5x{Xa]al) =iSxiX)^. . 
Eq. (52) can therefore be written 

(E + V)(° -;)X = VY, 

or 

ExpUcitly, these relations are 



(53) 



The column matrices and T-j- do not vanish because 
1$) is not an eigenstate of N^- or T±. Therefore (^°^ ^ 

and o)"'"t right eigenvectors of R with the eigen- 
values zero and ±/i. They belong to the tt and np spaces, 
respectively. For /x = 0, one has Ty = because (P) 
points into the y direction so that |$) is an eigenvector 
of Ty. The column matrices o)"'"=f ^^^^ merge into 
a single right eigenvector J^T^ of R with the eigen- 
value zero. 

For K = 0, the colomn matrices in the common right 

kernel of R and K found here must be those whose ex- 
istence was proved in Sect. IIIG. It is now seen that 
they persist for k 0. It follows from the discussion in 
Sect. Ill G that they span the kernel. The derivatives 
d'^R]iB/ dQi^dQ^ are for . = the elements of K. Since 



it has now been proved that K is positive semidefinite and 
its entire kernel stems from symmetries of -Rhb, it fol- 
lows that 1$) minimizes Rbb locally. For a valence space 
consisting of a single j-shell, and an energy functional 
without the symmetry force and including the complete 
expectation value of the pairing force, Camiz, Covello, 
and .lean [48] prove that a state with annihilators of the 
form (13) minimizes this functional globally on a cer- 
tain class of Bogolyubov quasinuclcon vacua with contant 
(iVV). Using their method, one can easily prove that |$) 
minimizes E^m globally on the same class of Bogolyubov 
quasinuclcon vacua. 

I promised in Sect. IIIC to prove that a quasinu- 
clcon vacuum annihilated by an operator and obeying 
Eq. (11) cannot minimize i?HB locally. Assume, to the 
contrary, that |$) is annihilated by an operator at. In 
order to satisfy Eq. (11), the number of such operators 
must be even for each kind of nuclcon, so for one kind of 
nuclcon there must be at least two such operators, 

and a^,^, say. These operators replace a^^ and a^,^ in 
Eq. (50). Matrices E and V can be defined as the Hessian 
matrices of (i?o) and i?HB — {Rq) with respect to Cik for 
Cfej = 0. In particular, 

= — Eqr — Eg'r < , 

while V is still given by an expression in the form of 
Eq. (40) with real P^, P^, T_|_, and T^. Now, let a column 
matrix x be defined by 

Xq'rgr = -X(qT)-i (9't)-i = 1 > X(.k = Otherwise. 

Because Tz is Hermitian, o)"'"« ~'^z' 

T^x = Tjx = . 

The contribution to x^ ^ J Vx from the symmetry force 

therefore vanishes. The contribution from the pairing 
force is nonpositive. With K given by Eq. (30), then 

x'^Kx < , 

so K is not positive semidefinite and |$) does not mini- 
mize i?HB locally. 

The eigenvalues zero of R in the tt spaces can be un- 
derstood to result from the ground state of R being an 
eigenstate of Nr- The similar applies to the eigenvalues 
zero in the np space for /z = 0. For ^ > 0, the eigenvalues 
±/i can be understood to originate in the M-j' = T ~ 1 
isobaric analog of the ground state, which has the excita- 
tion energy /i as an eigenstate of R. This interpretation 
differs from that of Ginoccio and Wesener [11], who, in 
the quasiboson picture, find that the vibrational quan- 
tum in question increases T by one unit. They infer this 
from the commutation relation, which they do not prove, 
written in the line after their Eq. (88b). 
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In the oscillator interpretation of the expression (39), 
the term i/i arises from the quantal fluctuations of the 
variables T±, which are in the approximation [T_|_,T_] = 
2Tz ~ 2T proportional to conjugate coordinates. In an 
eigenstate of and T^, these fluctuations are deter- 
mined by the isospin algebra. Assume, in accordance 
with the isorotational picture, that intrinsic classical vari- 
ables, invariant under isorotation, could be defined so 
that the Hamiltonian is a function of |T| and these in- 
trinsic variables, and let Ecoi{\T\), where "col" stands for 
"collective" , be the energy of the intrinsic equilibrium for 
given |T|, where T is a "classical" isospin. The term ^/i 
in Eq. (39) then correponds to the following estimate of 
the energy due to the quantal fluctuations of T. 



i^coi(Vl^(T+l)) - Scoi(T) 
AE 



A^coi(T) 
AT 



1 ^ 

2 AT 



where the differences are taken between the discrete al- 
lowed values of T and the T dependence of the total zero 

point energy is neglected in the second approximation. 
A related interpretation of the a term involving the an- 
gular velocity which emerges analogously in models with 
rotational invariance is discussed by Marshalek [44] in a 
boson expansion picture. 



Since T, 



2 



NJ and 



it follows from Eq. (32) that any other right eigenvector 
of R than those mentioned so far is orthogonal to the 
components of T. As seen from Eqs. (28) and (40), it 
then belongs to the kernel of the symmetry-force term 
in Eq. (40). Therefore, ±/x are the only eigenvalues of R 
which depend on k and the right eigenvectors of R are 
independent of k. Since, as seen from Eqs. (14), (15), 
(16), and (17), and are independent of k. it hence 
follows that it is sufficient to determine the eigenvalues of 
R for K = 0, in which case the nn and pp spaces separate. 
The symmetry force is then taken into account by adding 
the term kT in Eq. (18). 

This is a remarkable result. Its interpretation in the 
isorotational picture is that the symmetry force docs not 
influence the intrinsic excitations. It only contributes to 
the collective energy. This is consistent with the fact, fol- 
lowing from the remark after Eq. (6), that the symmetry 
force differs from ^kT^ only by a constant. 

Prom Eqs. (18), (19), (35), and (39) follow 

EupA = E^pA.K=o + ^kT , (54) 
E^Ehb+c + Erpa = E^=o + ^K(TiT + 1) - |Av) . 

(55) 



Using Eq. (45), = Jjj^, and the fact that 

o)"'"z belongs to the kernel of E -I- V, one can eas- 
ily prove that M' defined by Eq. (49) is independent of k 
in the nn+pp space. Thus, its elements between the nn 
and pp spaces vanish. 

I. Case = 

Most of the preceding discussion applies for A,- = 0. 

The Fermi level can be placed, in this case, anywhere 
between such a pair of single-nucleon levels that Nr/2 
levels Cq satisfy < Ar- The pairing energy -Epair van- 
ishes according to Eq. (21), and one gets from Eqs. (19) 
and (20) 



E, 



HB 



2E E + 



(56) 



The column matrices vanish, so R has no eigenvalue 
zero in the nn+pp space. The brackets (T_)_^^ also 

vanish, and (T_)^_i__i^ is equal to \/2 if Ap < < An 

and zero otherwise. For T = 0, the column ma- 
trix T^ then vanishes because no level satisfies 
Ap < Eg < An, so R has neither in the np space an eigen- 
value zero. For T > 0, consider the part of the np 

space spanned by the operators [a^a^]^__|_ and [aa]_^_|_ 
with Ap < < An. In this subspace, (^^^ o)"'"- 

the elements \/2 in its upper part and zero in its 
lower part. From (Pj)_^^ = (Pz) = follows 

R = for K = 0, where the diagonal elements of 

E are Eqn + Eqp = (An - eq) + [eq - Ap) = An - Ap. 
Thus o)"'"- for K = an eigenvector of R with the 
eigenvalue An — Ap = /U. This is the only positive eigen- 
value of R which depends on k. For k > 0, it becomes 
An - Ap -I- kT = /X. 

For G = K = 0, one has ^^rpa = 0. Generally for 
G = 0, Eq. (54) then gives 



Erpa = ^kT . 



(57) 



Eqs. (35), (55), (56), and (57) imply Eq. (6), so the RPA 
is exact in this case. For G > 0, the pairing force gives 
an additional contribution to i^RPA. 



J. Case of particle-hole symmetry- 
Special phemonema occur when the single-nucleon 
spectrum is symmetric about some energy and the va- 
lence space is halfway filled. For brevity, I call this 
symmetry of the single-nucleon spectrum a "particle- 
hole symmetry". Without loss of generality, the spec- 
trum can be assumed, in this case, to be centered at 
zero and labeled so that = —e-q > for gr > 0. Then, 
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A = 0, At: 



-An, A, 



An , E„ 



Vqp = U(_q)„, and M' defined by Eq. (49) has equal eigen- 
values in the nn and pp spaces. 

In the np space, the eigenvalues of M' turn out to have 
the multiplicity 2 with the exception of the eigenvalue /x^ 
found in Sect. IIIH and an eigenvalue 4(A^ + A^), both 
of which are for k = less than all {Eg^ + Egp)"^ and 
less than all other eigenvalues of M'. To see how this 
happens, first consider the case k = 0. Let arr' i ^tt' be 
divided into blocks corresponding to an ordering of the 
subscripts q by their signs with g > first and a common 
order of \q\ for each sign. Eqs. (47) and (48) then give 

i (i -i) (^^P + bnp)(anp - b„p) (I 



e 

e-2Gp_pl 







1/2 T 1/2 

2Ge ' P+P+e ' 



- 2Ge 



1/2 T 1/2 

' p pe ' 



where the similarity is effected by the matrix 



V ei/V' 



-1/2 Q 

1/2] ■ The matrices e and p± are the upper left 

and upper parts of Cnp and Pnp±- From Eqs. (14), (15), 
and (48), one gets by some algebra 

{Eq^ + i^9p)Pg„p± 
= l(^ + ^] ((^.n + Eq,r - A±) (58) 



E E 



with 



A+ = 4A2, A_=4{Xl + Al) 



(59) 



It follows from Eq. (16) that the column matrix s with 
the dimension of Pnp± and the elements 



(60) 



is a unit vector. Prom Eq. (14) one can derive 
'qnVqp for Q > Q. It thcu foUows from Eq. (48) 
that the elements of p± are positive. It follows from 
Eq. (58) that A± are then less than all {Eg^ + Egp)^, so 
the diagonal matrices 

f± = e^-A± (61) 
are positive definite and the positive eigenvalues of 

g = (l-ssV(l-ss^) 

are greater than A± . 

From Eqs. (58), (60), and (61), one gets 



A. 



It is seen that s is an eigenvector of h± with the eigenvalue 
A±. Furthermore, if > is an eigenvalue of g and x 
the corresponding eigenvector, then 



A± 



T 2 

ss e x 



is an eigenvector of h± with this eigenvalue. The ma- 
trices h-t thus have these eigenvalues in common. The 
eigenvalue A+ = (2An)'^ = /i^ of M' is just that discussed 
in Sect. IIIH. It is the only eigenvalue of M' which de- 
pends on K. For k > 0, it becomes (2An + kT)"^ = ji^. 
The eigenvalues of h± are seen to satisfy 



T 

P±;2 



2e 



1 

G 



(62) 



Thus, M' has in the np space one twofold degenerate 
eigenvalue in each interval between consecutive squares 



of 



Eqp, q>0. 



For T = 0, the eigenvalues of M' in a rr space and 
their multiplicities are the same as in the np space. The 
eigenvalues with the multiplicity one are zero and 4A^ 
in this case. It follows that, more generally, if the energy 
levels of the nucleons with isospin mt are symmetrically 
distributed about A^, then the eigenvalues of M' in the 
TT space have the multiplicity two with the exception 
of eigenvalues zero and 4A^, which have the multiplicity 
one. This was proved previously for k = by Hogaasen- 
Feldman [49] and Bes and Broglia [50]. 



IV. CALCULATIONS 
A. Equidistant single-nucleon levels 

An insight into the general behavior of the model 
is obtained by studying the idealized case of infinitely 
many equidistant single-nucleon levels. This is the topic 
of the present section. More precisely, I assume that 
the system has d valence single-nucleon levels e,, spaced 
by a constant D and symmetrically distributed about 
zero, and that = 2d. The discussion in Sect. IIIJ 
then applies, and I denote by A the value of An for 
T = 0. This system is considered in the limit d — » oo 
with D, A, and k kept fixed in the limit. Appropri- 
ate values of D, A, and k can be derived from empir- 
ical formulas in the literature. Thus, D = Aj (pa/TT^) 
with a = 0.176A(l - l.OA''^/^) MeV"^ according to 
Kataria, Ramamurthy, and Kapoor [51]. An empiri- 
cal formula of Bohr and Mottelson [3] for the odd-even 
mass difference gives A = 12 MeV. The con- 
stant (D -I- k)/2 will be seen to be close to the coeffi- 
cient of T{T + 1) in a mass formula. According to an 
empirical mass formula by Duflo and Zucker [52], then 
D + K = 2(134.4^-1 - 203.6A-4/3) MeV. Table I shows, 
for some mass numbers A, the values of D, A/D, and 
k/D given by these expressions. 

Since the dependence of i? on k is trivial according 
to Eq. (55), I assume for now k = 0. In units of D, 
the energy E is then a function of A/D and T. Some 
quantities mentioned in the following are infinite in the 
limit d^ oo. When such quantities are said to obey some 
relation, this is meant to be be arbitrarily accurately true 
when d is sufficiently large. I have checked that these 
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TABLE I: Parameters and results of the calculations with 
infinitely many equidistant single-nucleon levels. 



A 


D (MeV) 


A/D 


k/D 


a 


b 


X 


48 


1.1 


1.6 


2.0 


.873 


4.08 


1.010 


100 


.5 


2.5 


2.8 


.859 


5.74 


1.008 


240 


,2 


1.2 


:!.() 


.812 


8.71 


1.005 



statements hold very accurately in calculations for d = 
500 and A/D = 2.5. 

The Hartree-Bogolyubov energy Ehb is according to 
Eq. (19) the sum of Eq and impair, both of which are 
infinite in the limit d oo. In this limit, the only change 
in the solution of Eq. (16) and (17) when Nr changes 
by two units is a change of A,- by the level distance D. 
Therefore, = A independently of T, so impair does not 
contribute to the symmetry energy. On the other hand, 
Eq docs so. This is because the isospin T is produced by 
the promotion of T nucleons from proton levels with the 
average energy Ap/2 to neutron levels with the average 
energy An/2. The result is an increase 

P F — An ~ Ap ^ 2 



where 

^^=A„ = -Ap=T^ = ii)T 

has been used. 

Eqs. (39), (59), and (63) give 



(63) 



ErpA = Enn + En 



C,DT) 



2 .2 
+ A , 



(64) 



Enn — ^ ^ ^nn,n ^ ^ ^ Eg^ , 
n q 

Enp — ^ ^ ^np,n ^ ^ Eqn , 



where Wnn,n are the positive eigenvalues of R in the nn 
or pp space and C0np,n the twofold degenerate, positive 
eigenvalues of R in the np space. By adding the third 
term in the expression (64) to Eq, one gets a contribution 
to the symmetry energy proportional to T(T + 1), 

Eo - Eo,T=o + \DT = \DT{T + 1) . 

The energies £Jnn and E^p are infinite in the limit ^ 00 . 
Because the quasinucleon energies Eq^ and the elements 
of R in the nn space arc just relabeled when is changed 
by two units, E^n is independent of T. On the other 
hand, E^p gives a nonzero contribution to the symmetry 
energy. This is demonstrated in Fig. 1 by a calculation 
for A/D = 2.5 and d = 1000. A calculation for d = 500 
gives a curve which cannot be distinguished from this 
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FIG. 1: The contribution of £np to the symmetry energy in 
units of D for A/D = 2.5 and d = 1000. 



one. The limit d ^ 00 is thus realized in practice for d 
of this order. 

The decrease of Enp with increasing T is under- 
stood from the fact noticed in connection with Eq. (62) 
that the frequencies tOnp^n lie between consecutive two- 
quasinucleon energies Egn + Egp, q>0. The lowest 
two-quasinucleon energy increases with T, whereas 
when Eqn + Egp is sufficiently large, Eqn + Eqp « 
(eg — An) + {eg — Ap) = 2eq, which is a constant. There- 
fore the difference between the sum of Wnp.n and the sum 
of Eqn + Eqp, Q > 0, that is, J2q Eqn, dccreases. The sum 
-E'np + v^Ai increases because all the positive roots u of 
Eq. (62), including \/A±, and all the two-quasinucleon 
energies arc bounded below by y/ A±, which increases. 
Therefore also E^rpa, which differs from the average of 
Enp + i/A-t only by the constant Enn, increases. 

The calculation of E^p for A/D = 2.5 and d = 1000 
was extended to larger T than shown in Fig. 1, up to 
T = 100. The entire set of results is well decribed by the 
expression 



E, 



np 



E, 



np,T=0 



-D(V(aT)2 + 62-6), 



a = .859 , b = 5.74 . 



In the view of Fig. 1, this gives a curve which can hardly 
be distinguished from the calculated one. Similar results 
are obtained for A/D = 1.6 and A/D = 4.2 with a and 
b given in Table I. 

Now collecting all the contributions discussed above 
and adding the last term in Eq. (55), one gets 



E~E' 



T=0 



i(D + K)T(T+l) 



-D ^{aTY + h^ 




d) -'^d 



(65) 
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For the paramtcrs in Table I, the term subtracted from 
^(-D + k)T{T + 1) in this expression is positive for T > 0. 
It amounts to 11% of ^(-D + k)T for ^ = 48 and T = 8 
C**S), 10% for A = 100 and T = 14 (i""Kr), and 9% for 
A = 240 and T = 28 (^^Ou). The hnear term \{D + k)T 
is thus the dominant correction to the quadratic term 
i(D + k)T'^. The sign of the additional term in Eq. (65) 
is, however, consistent with the experience, mentioned in 
Sect. II, that in global fits to the empirical masses with 
formulas which include a term proportional to T(T + x), 
the constant x tends to be somewhat less than one. For 
T w 0, Eq. (65) gives 

X = 1^1 - j • (66) 

Table I shows x calculated from Eq. (66) for the param- 
eters in the table. It is seen that x « 1.01 in these cases. 

As discussed in Sect. Ill I, the RPA gives the exact 
energy (6) for G = 0. The symmetry energy derived 
from this expression in the present case is 

E - Et=o ^^{{D + k)T^ + kT) , (67) 

without the term \DT. The spontaneous breaking of the 
isobaric invariance by the pairing force is thus required 
for this term to appear. 

B. Comparison with Skyrme force models 

Nuclear masses are often compared with Hartree- 

Fock (HF), Hartree-Fock-BCS (HFBCS), or Hartree- 
Fock-Bogolyubov (HFB) calculations with phenomeno- 
logical energy functionals based on Skyrme forces. The 
HFB method is described in detail in a recent article 
by Chamel, Goriely, and Pearson [53], where also refer- 
ences to earlier work in this line of research are found. 
Applying the HF, HFBCS or HFB scheme in an approx- 
imate way to the present Hamiltonian sheds a light on 
the origin of certain phenomena observed in the Skyrme 
force calculations. In the formalism of Ref. [53] , a pair- 
ing and a particle-hole part of the two-nucleon interaction 
are treated differently, as in the present theory. For the 
present pairing interaction and (P^) ~ 0, there is then 
no difference between the HFBCS and HFB schemes. 

Apart from the different interactions, the only differ- 
ence between the formalism of Sect. IIIC and the HFB 
formalism of Ref. [53] is that the particle-hole matrix el- 
ement is antisymmetrized in the latter. In the formalism 
of Sect. Ill C, this amounts to replacing i?HB with 

-Rhb - \Ktxt ■ ptp, {j\p\k} = {a\a^) . (68) 




5 10 15 
T 



FIG. 2: The quantity 9 - 6»t=o as a function of T for d = 1000 
and A/D = 2.5. 



The resulting Routhian is stationary at the quasinucleon 
vacuum |$) given by Eqs. (7), (13), (14), (15), (16), and 
(17) when is replaced with €q — jK{Vg^ + 2v^^,), t' ^ r. 
I neglect this modification of the self-consistent single- 
nucleon energy and calculate the second term in Eq. (68), 
the "Fock term", with Ugr and Vgr given by Eqs. (14) 
and (15). The essential factor is 

trt • ptp = 5](i(<n + <p) + Hn^lp) 
9 

3 1 _ ( 1 ( 2 2,2 2\,2 2,2 2 \ 

(69) 

The first term in this expression gives the second term in 
the outer bracket in Eq. (35) and the second term a part 
of the first of the diagrams (36). In the case of infinitely 
many equidistant singlc-nuclcon levels, the only part of 
the expression (69) which depends on T is 

S = Y.i<y^V+«r)- (70) 

9 

The quantity 9 — 9t=o is shown in Fig. 2 as a function 
of T for d = 1000 and A/D = 2.5. It is seen that the 

asymtotic slope is one. This is easily understood from 
Eq. (70). Indeed, for T ^ 1 the first term in the bracket 
is practically zero, while the second term is close to one 
for Ap < Eg < An and close to zero otherwise. Hence, 
for T ^ 1, the contribution of the Fock term to the 
symmetry energy deviates from the linear term in the 
expression (67) only by a constant. For T « 0, it is, how- 
ever, quadratic in T because is by Eqs. (14), (15), (63), 
and (70) analytic and even as a function of T. Therefore, 
the HFB scheme does not produce any cusp at iV = Z in 
the curve of masses along an isobaric chain. The study 
by Satula and Wyss [22] shows that this is true also in 
HFBCS calculations with Skyrme forces, and it is true 
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in the HFB calculations discussed in Ref. [53], where, 
in order to reproduce the empirical masses, a so-called 
Wigner correction is added to the HFB energy. 

HF calculations with Skyrme forces correpond to the 
case G = of the present theory. 1 this case, the 
Fock term is equal to ^k{T — jA^), so the exact en- 
ergy (6) is recovered. With infinitely many equidistant 
single-nucleon levels, the symmetry energy is given by 
Eq. (67) and thus includes the linear term ■^kT. Satula 
and Wyss [22] find that such a term appears also in HF 
calculations with Skyrme forces. As in the present the- 
ory, there is in these calculations no term corresponding 
to ^DT. The present schematic model thus seems rep- 
resentative of the Skyrme force models as to the basic 
mechanisms responsible for the absense or presense of 
such linear terms of various origins. 



C. Deformed Woods-Saxon single-nucleon levels 

Moving now from the idealized case of a uniform single- 
nucleon spectrum to assumptions closer to the reality, I 
discuss in this section calculations with single-nucleon 
levels derived from a deformed Woods-Saxon potential. 
I have chosen for this study isobaric chains whoso T = 
mcmbcirs have low-lying 2"*" levels and can thus be sup- 
posed to have appreciable quadrupole deformations. For 
each isobaric chain, the single-nucleon energies are calcu- 
lated for a deformation pertaining to the T = nucleus. 
The definition of the deformed Woods-Saxon potential 
follows Dudek et al. [54]. The average nucleon mass and 
the average of the "universal" neutron and proton param- 
eters of Dudek et al. [55] are employed, and no Coulomb 
potential is included. The halfdepth surface of the spin- 
independent part of the potential is assumed to have a 
pure, prolate quadrupole shape with the deformation /3 
derived by the relations of Raman, Nestor, and Tikka- 
nen [56] from B{E2,0~^ 2+), when it is known, and 
else from the 2+ excitation energy. All bound single- 
nucleon states arc included in the valence space. 

The pairing force contant G is determined for each 
isobaric chain by demanding An = Ap = 12A~^/^ MeV 
for T = 0. Most other nuclei in the isobaric chain then 
acquire similar values of A^, but A^ may vanish if 
reaches a gap in the single-nucleon spectrum. The calcu- 
lated symmetry energy turns out to be minor sensitive to 
the precise value of G. The constant k is chosen so as to 
best fit the data. A summary of the resulting parameters 
is given in the first, third, and fourth row of Table II. A 
somewhat disturbing feature of these parameters is the 
slightly irregular A dependence of the optimal k. This 
irregularity may be related to the crude treatment of the 
shape degrees of freedom. 

For comparison with the calculated results, I have de- 
rived an empirical symmetry energy from the masses 
compiled by Audi, Wapstra, and Thibault in their 2003 
Atomic Mass Evaluation [57]. It is calculated from the 
binding energy for Mt > minus the electrostatic term 



TABLE 11: Parameters and results of the calculation with 
Woods-Saxon single-nucleon levels. 

A (3 G (MeV) Interval of (MeV) k (McV) 



48 .337 .450 1.62-1.85 1.4 

56 .4 0-1.94 1.1 

68 .234 .286 0-1.65 1.5 

80 .342 .262 1.08-1.69 1.1 

100 .2 0-1.55 .8 



extracted by Kirson [29] from the differences of the bind- 
ing energies of mirror nuclei, 

-.716^2-^.993^4/3 

Am 

The calculated and empirical symmetry energies are 
shown in Fig. 3. In view of the schematic character of 
the model, where, in particular, any variation of shape 
degrees of freedom is absent, the agreement is suprisingly 
good. Also shown in the figure is the composition of the 
calculated symmetry energy. The term ^kT^ typically 
makes up about half of it for large T. The contribution 
of the pairing energy impair is plus or minus a few MeV. It 
arises from the variation of Ar along the isobaric chain. 
This variation is due, in turn, to the variation with At of 
the local single-nucleon level density. The irregularity of 
the single-nucleon spectrum is also responsible for a fairly 
irregular behavior of the contribution from Eq as com- 
pared to its quadratic dependence on T in the idealized 
case. 

The part of the contribution from E^pj^ originating in 
the nn+pp space is essentially zero, as in the idealized 
case. The part originating in the np space is dominated 
by the single eigenvalue ^/z of R, which increases essen- 
tially linearly with T. Its deviation from ^/i appears to 
be of higher than linear order in T and always negative, 
just as in the idealized case, and it amounts to about 
— 1.5 MeV for the largest T's for which the binding energy 
is measured. The ratio (-Erpa — -Brpa,t=o)/(-E — Et=o) 
is close to T/(T2 -hT) = 1/3 for T = 2 in the 
three cases. This is understood from the follow- 
ing facts. (1) i?HB — -£'HB,(T^)=o oc (Tj;)^ for small 
(T^) because E'hb is analytic and even as a fun- 
tion of {Tz). (2) Srpa — i?RPAT=o is dominated by 
y=\dEnB/d{T,). (3) (T,) =r.' 

D. Spherical Woods-SeLxon single-nucleon levels: 
yet another mechanism 

For ^ = 56 and A = 100, the T = nucleus is dou- 
bly magic. Therefore, I have made calciilations for these 
mass numbers with a spherical Woods-Saxon potential. 
I have chosen G in this case so that A,- 12^"^/^ MeV 
for T > 0. Then, A^ = for T = 0. The adopted param- 
eters are given in the second and fifth row of Table II, and 
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FIG. 3: Calculated (solid line) and empirical (crosses) 
symmetry energy for some values of the mass number A. 
Shown with different signatures are also the following com- 
ponents of the calculated symmetry energy: Eq — So,t=o 
(dash-dot), -Epair — iJpair,T=o (dash-dot-dot) , |kT^ (dash- 
dot-dot-dot), i?RPA,nn+pp - -BR,PA,nn+pp.T=() (short dash), 

^RPA.np — iSRPA,np,T=o (loug dash), where iSRPA,nn+pp and 
£^RPA,np are the contributions to the expression (39) from the 
nn+pp and np spaces. The dotted curve is i/it. 
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FIG. 4: See the caption to Fig. 3. 



the results of the calculations are shown in Fig. 4. The 

valence space consists of the levels below the Nr = 50 
shell closure for A = 56, and below the Nt = 82 shell 
closure for A = 100. 

The contributions to the symmetry energy from i kT^ 

and -Erpa are similar to those of the deformed case. But, 
due to the increase of when A,- moves into a shell 
above or below the magic gap, Spair gives in this case 
a significant negative contribution. Because the largest 
decrease of -Bpair takes place for low T, E'pair then con- 
tributes to the nonlinear part of the symmetry energy. 
On the other hand, Eq is practically linear. This is be- 
cause the increase of Eq results from the promotion of 
nucleons from proton states below the shell gap to neu- 
tron states above the gap. The slope of Eq is thus roughly 
equal to the shell gap. The slope of Eq + £^pair at low T 
is similar to that of -Erpa for A = 56 and several times 
that of i^RPA for A = 100. In both cases, it thus gives a 
significant contribution to the linear part of the symme- 
try energy. This could be expected to be true in general 
in isobaric chains with a doubly magic T = nucleus. 
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V. EXPERIMENTAL SIGNATURE OF 
SUPERFLUID ISOROTATION 

An experimental signature of a nuclear rotation in 
space is the sequence of large and approximately equal 
reduced E2 transition probabilities between consecu- 
tive members of a rotational band produced by the 
quadrupole deformation of the intrinsic charge distribu- 
tion. The quadrupole deformation is measured by the 
mass or charge quadrupole tensor. In the case of super- 
fluid isorotation, the intrinsic deformation is measured 
by the isovector P. Since the consecutive members of an 
isorotational band are separated by two units of isospin, 
they are not connected directly by an isovector. However, 
as pointed out by Bohr and Mottelson [3] , the superfluid 
isorotational bands participate in a larger structure con- 
sisting of the ground states of all the doubly even nuclei 
and their isobaric analog states. The groimd states of the 
doubly even isotopes are connected by and the ground 
states of the doubly even isotoncs by Pp. The isovector 
component Pz connects such states with \Mt\ = T — 1 
isobaric analog states in doubly odd nuclei. The chains 
of superfluid isotopes or isotoncs in fact form the pair 
rotational bands discussed by Bohr [58]. Yoshida [59] 
shows that superfluidity enhances the ground state to 
ground state cross section of two-ncutron or two-proton 
transfer between doubly even nuclei by a factor about 
(2At/G)^, where r = n or p. Isovector one-neutron-one- 
proton transfer between a ground state of a doubly even 
nucleus and an isobaric analog state is then similarly en- 
hanced. 

The experimental signature of superfluid isorotation 
therefore coincides with that of pair rotation. The pic- 
ture of a superfluid isorotation implies that the enhance- 
ment factors of two-nucleon transfer involving doubly 
even ground states or their isobaric analogs remains ap- 
proximately constant all the way down to T = 0. Near 
closed shells, the pair rotational bands may develop into 
the pair vibrational bands dicussed by Bohr [58], which, 
as well, have enhanced two-nucleon transfer cross sec- 
tions. If, on the other hand, doubly even nuclei with 
T = would have a structure radically different from that 
of doubly even nuclei with T > 0, a major deviation from 
such a smooth behavior would be seen. No such effect 
seems to be indicated by the enhancement factors com- 
piled by Bes et al. [60]. 



VI. SUMMARY 

A Hamiltonian with a single-nucleon term and a two- 
nucleon interaction was investigated. The single-nucleon 
Hamiltonian has fourfold degenerate eigenvalues corre- 
ponding to time reversed pairs of neutron and proton 
states. The two-nucleon interaction has a pairing and 
a particle- hole part. The paring part is the isobari- 
cally invariant isovector pairing force and the particle- 
hole part an interaction of isospins, which I call the 



symmetry force. A Routhian was contracted by sub- 
tracting from the Hartrce-Bogolyubov energy functional 
Lagrangian multiplier terms propotional to the expecta- 
tion values of the number of valence nucleons and the 
z component of the isospin, and it was shown that this 
Routhian is locally minimized by a product of neutron 
and proton Bardeen-Cooper-Schrieffer [23] states. 

This quasinucleon vacuum and a single- quasinucleon 
Routhian operator derived from the Hartree-Bogolyubov 
Routhian was taken as the starting point for a calcu- 
lation of the ground state energy as a function of the 
number of valence particles and the isospin quantum 
number T in the Random Phase Approximation (RPA). 
The correction to the Hartree-Bogolyubov energy is the 
sum of a term which does not depend on T and a term 
Ekpa equal to half the sum of the poles of the RPA 
two-quasinucleon propagator minus the sum of the two- 
quasinucleon energies. The poles of the two-quasinucleon 
propagator which are different from two-quasinucleon en- 
ergies can be determined separately in a two-ncutron, 
a two-proton, and a neutron-proton quasiparticle space. 
In each of these spaces, there is a Nambu-Goldstone pole 
due to the global gauge invariancc and isobaric invariance 
of the Hamiltonian. The two-neutron and two-proton 
Nambu-Goldstone poles have the frequency zero, while 
the neutron-proton Nambu-Goldstone pole is equal to the 
Lagrangian multiplier of the z component of the isospin. 
The term in i?RPA resulting from this pole was inter- 
preted in a picture of a collective rotation in isospace to 
be due to the quantal fluctuation of the isospin. The pole 
in question is the only one which depends on the strength 
K of the symmetry force. The only contribution of the 
symmetry force to the symmetry energy is therefore a 
term ^kT{T+1). 

If the single-nucleon spectrum is symmetric about a 
certain energy and the valence space halfway filled, the 
neutron-proton poles are twofold degenerate except for 
the Nambu-Goldstone pole and one more pole, which 
has an analytic expression. Related results are known 
from the literature to pertain to the neutron-neutron and 
proton-proton spaces. 

In an idealized case of infinitely many equidistant 
single-nucleon levels, the neutron and proton systems 
have a common pair gap A which does not depend on 
T. Therefore, the pairing force does not contribute to 
the symmetry energy. Neither do the two-neutron and 
two-proton parts of -Erpa- For k — 0, the single-nucleon 
term in the Hamiltonian and the neutron-proton Nambu- 
Goldstone pole give together a contribution equal to 
^DT{T -1-1), where D is the single-nucleon level spac- 
ing. The second nondegenerate neutron-proton pole gives 

a contribution equal to (^DT)^ + — A. The re- 
mainder of the contribution of the neutron-proton part of 
Erpa was calculated numerically. In a very good approx- 
imation, it is —d(^-^ (aT)"^ + 6^ — where a and b are 

fimctions of A/D. For realistic parameters, the sum of 
these two terms is negative and amounts to about —10% 
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of the linear term ^{D + k)T for the largest T of ob- 
served nuclei. For T w 0, they give a contribution to the 
symmetry energy quadratic in T which makes the sym- 
metry energy proportional to T{T-\-x) with a; « 1.01. In 
the absence of the pairing force, the RPA gives the ex- 
act symmetry energy, which is i ((D -|- k.)T'^ + kT) . This 
expression does not have the linear term \DT, which 
thus appears only when the isobaric invariance is spon- 
taneously broken by the pairing force. 

If the matrix element of the symmetry force is anti- 
symmetrized and the contribution of the exchange term 
to the self-consistent singlc-nucleon energy and the RPA 
correlations are neglected, the total contribution of the 
symmetry force to the symmetry energy is asymptoti- 
cally for large T equal to ^kT{T + 1) plus a constant. 
For T « 0, it is quadratic in T. Therefore the curves 
of masses along isobaric chains get no cusps at T = 0. 
This corresponds to observations reported from Hartree- 
Fock-Bogolyubov calculations with Skyrme forces. In the 
absence of the pairing force, the exact symmetry energy 
without the linear term ^DT is recovered also in this ap- 
proximation. This correponds to observations reported 
from Skyrme force Hartrce-Fock calculations. 

Calculations with Woods-Saxon single-nucleon levels 
give results in surprisingly good agreement with the em- 
pirical variation of the binding energy of doubly even 
nuclei along isobaric chains. In the case of a deformed 
Woods-Saxon potential, the behavior of the individual 
components of the calculated symmetry energy is similar 
to their behavor in the idealized case. In calculations for 
A = 56 and A = 100 with spherical Woods-Saxon levels, 
the promotion of nucleons across the N — Z = 28 and 
N = Z = 50 gaps in the single-nucleon spectrum gives a 
large linear contribution. Due to the onset of superfluid- 
ity when the neutron and proton Fermi levels move into 
the shells around these gaps, the pairing force gives an- 
other large contibution in this case. Together, these two 
contributions give a linear term which is comparable to 
or larger than that of i^RPA- 

In a picture of a collective rotation in isospace, the 
intrinsic deformation is measured by the pair annihila- 
tion isovector. The isorotation can therefore be char- 
acterized as superfluid. The pair annihilation isovector 
does not connect directly consecutive members of a su- 
perfluid isorotational band, which differ by two units of 
isospin. However, the superfluid isorotational bands par- 
ticipate in a larger structure which includes the pair ro- 
tational and pair vibrational bands discussed in the lit- 
erature. Within these bands, the cross sections for two- 
nucleon transfer are enhanced by the superfluid correla- 
tions. The picture of a superfluid isorotation implies that 
the enhancement factors of isovector two-nucleon trans- 
fer should develop smoothly down to T = 0. This seems 
to be consistent with the empirical evidence. 
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APPENDIX: WIGNER ARGUMENT FOR THE 
ISOBARIC SU{2) 

Consider a system of valence nucleons and a two- 
nucleon interaction X^i;, where the summation is over 
pairs of nucleons, and let Pg and Pa project to the spaces 
of two-nucleon states symmetric and antisymmetric in 
position and spin. In any state of the system, 

<E«) = (E(^s + PaH 

The total antisymmetry implies Pg = ^ — ti • t2, where 
ti and t2 are the isospins of the two nucleons. Hence 

_ iv(iv + 2) 
^ '~ 8 2 ' 

EPa = E(i-Ps) = '^-^^;"'^ + ^. 

The assumption in the case of the isobaric SU{2) sym- 
metry corresponding to Wigner's in the case of the 
SU{4:) symmetry is that the ground state values of 
(E Psv)/{J: Ps) and )/(E -Pa) are smooth fimc- 

tions of Ay and T. If, in particular, they do not depend 
on T, the symmetry energy is proportional to T(T +1). 

The assumption of T-independent Psv) / {J2 Ps) 
and (E-P*at')/(E-fa) is easily seen to be invalid in 
the case of the isobarically invariant isovector pair- 
ing force acting in a single j-shell. In this case, 
V = —{2j -h l)Gi-bi where Pq projects to two-nucleon an- 
gular momentum zero. Prom an expression for the eigen- 
values of (2j -I- 1) E ^0 derived by Edmonds and Flow- 
ers [61], it follows that the lowest eigenvalue oi^v for 
fixed even A^ and Av/2 -|- T is 

G{TiT+l)-^:^^l±^^). (A.2) 

Since a state of two nucleons from the same j- 
shell with angular momentum zero is antisymmetric 

in position and spin, and v is negative scmidefinite, 
P^v) = and {Y, P^v) < unless {J2v) = 0. Tak- 

i^S (E Pa-^)/{'l2 Pa.) to be a negative constant, one gets 
from Eq. (A.l) 

with a negative constant of proportionality. This expres- 
sion obviously conflicts with the expression (A.2). It even 
differs from the latter by the sign of its contribution to 
the symmetry energy. 
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